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Abstract. The model of competition between densities of two different species, called 
predator and prey, is studied on a one dimensional periodic lattice, where each site 
can be in one of the four states say, empty, or occupied by a single predator, or 
occupied by a single prey, or by both. Along with the pairwise death of predators and 
growth of preys, we introduce an interaction where the predators can eat one of the 
neighboring prey and reproduce a new predator there instantly. The model shows a 
non-equilibrium phase transition into a unusual absorbing state where predators are 
absent and the lattice is fully occupied by preys. The critical exponents of the system 
are found to be different from that of the Directed Percolation universality class and 
they are robust against addition of explicit diffusion. 



PACS numbers: 64.60.ah, 64.60.-i, 64.60.Dc, 89.75.-k 
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Absorbing configurations do not have any outgoing rates [T]. Once reached 
there, the system can not escape from these configurations. Presence of absorbing 
configurations in a phase space raises a possibility that the concerned system may 
undergo a non-equilibrium phase transition into absorbing states. The critical behavior 
of these absorbing state phase transitions (APT)s [2] depends on the conservation in 
dynamics and the symmetry between absorbing states. It has been conjectured [3] 
that in absence of any special symmetry the APT belong to the directed percolation 
(DP) universality class as long as the system has a single absorbing state. Additional 
symmetries, like particle- hole symmetry [J], conservation of parity [5], and symmetry 
between different absorbing states [6] lead to different universalities. Spreading process 
with spatially quenched randomness [7] or with long-term memory [8] are known to 
destroy the critical behavior completely, whereas the long-range interaction leads to 
continuous variation [9] of critical exponents. Presence of infinitely many absorbing 
states may [10] or may not [11] belong to DP universality class. Again a different 
critical behavior is observed when the activity field does not have any special symmetry, 
but it is coupled to a conserved density [12]. Recent studies have indicated that DP- 
critical behavior is possible, even in presence of an additional conserved field [13] . It is 
not quite clear, what microscopic ingredients can make an APT belong to the DP class. 

The models of directed percolation has been extended to more than one species |14j . 
Along with the simple DP behavior, a line of first order transition [15] has been observed 
in 1 + 1 dimension when two species compete for survival. Janssen [16] studied coupled 
DP processes with bilinear and bidirectional interspecies couplings in the framework of 
bosonic field theory, where no other critical phenomena were found other than the DP. 
Hierarchy of unidirectionally coupled DP processes with many species show multicritical 
behavior [17]. Coupled percolation processes have been also studied [18], where the 
absorbing phase become unstable with respect to an arbitrarily small branching rate 
even in one dimension. 

Predator-prey cellular automaton models [19] in two dimension show DP 
universality class. Coupled directed percolation (DP) processes with more than two 
species of particles (in one dimension) with different kind of interspecies coupling have 
shown DP-type [20] transitions. Lotka-Volterra like models in one dimension always 
show coexistence [21], either in form of well mixed states or as irregular bursts of the 
predator and prey population. A four state predator prey model [22] in one dimension 
with a restriction that a site can have at best one particle of each kind, shows an APT 
to an absorbing (extinct) state which belongs to DP-class. 

In this article we study a model of two species, say A (prey) and B (predator), 
on a (1 + l)-dimensional lattice. Each lattice site is either vacant or occupied by at 
best one particle of each kind. The preys grow independently as A0 — > AA and the 
predators die as BB — > 00, whereas they interact through a process BA — > BB, where 
birth of a new predator occurs instantly along with the death of the prey. The system 
show a line of continuous absorbing state transition different from DP as the rates of 
these processes are tuned. Unlike other multispecies models, in the absorbing state 
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Figure 1. Schematic diagram of the 4-state predator-prey (4SPP) model 

both preys and predators do not separately extinct, rather predators extinct and preys 
proliferate to fill up the whole prey branch. The model is defined on a one dimensional 
periodic lattice with lattice sites labelled by i — 1, 2 . . . L. Each site i can be either 
vacant or occupied by a single particle A (prey), or occupied by a single particle B 
(predator) or by both particles (co-existing A and B), thus the model can be treated as 
4-state predator-prey (4SPP) model. More than one particle of any kind is not allowed. 
These hardcore restrictions on individual particles, where co-existence is allowed, can 
be realized alternatively by considering two separate branches, one for A and the other 
for B particles. Particles living in one branch can not move to the other branch, as 
schematically shown in the Fig. [TJ Correspondingly, each site % is associated with four 
states; s l A — 0, 1 and s l B = 0, 1, where 1 (0) denotes the presence (absence) of a particle 
at site i. 

On a periodic lattice, these particles interact following a random sequential 
dynamics given below. The prey (A) can grow on their own branch with rate 
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p independent of the predators (5). Again, two neighboring predators (B) 
die simultaneously with rate q independent of A due to their own crowding or 
overpopulation. Here X in A-branch (5-branch) corresponds to an arbitrary state 
of A (B) particles i.e. presence or absence. These two species interact with rate r as 
follows; when a predator B at site i meets a prey A as the right neighbor with site 
i + 1 is not already occupied by a predator B, then the predator B eats the prey A and 
simultaneously reproduces another predator B. Note, that the dynamics is spatially 
asymmetric as both species grow only in one direction (here, right), and it does not 
include explicit diffusion of species. Effect on addition of symmetry and diffusion are 
briefly discussed towards the end. 

Let the steady state density of A and B particles be p s A and p s g respectively. Clearly 
in absence of predators (p s ) the prey density p s A = 1 as the preys do not have an 
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independent death process. Again, note that the dynamics do not allow A particles to 
be eaten when they co-exist with B on the same site. This indicates that these As can 
only die after their coexisting Bs die, which can occur with rate q if there is a B particle 
present or created at the immediate left neighbor. Thus, for reasonably small death rate 
of predators, p s A is expected to have a small value (p 8 < 1) when p s B ~ and then it 
increases along with p B ; so, the prey density p s A can never vanish. The predator density 
p s B can, however, become zero by repetitive death process. The isolated Bs wait until 
the prey invades their neighboring site and then they subsequently eat and reproduce 
with rate r and die with rate q. So, along with the coexisting phase [j] where both p s A 
and p s B are non-zero, we have another phase where p 8 = 1 and p s B = 0. Clearly, the 
later phase (p 8 = l,p s B = 0) is absorbing as once all the predators die, even then the 
single surviving prey can lead to proliferation of the prey population in the whole lattice. 
Of course, as argued earlier, p s A can not vanish, and the other possible absorbing state 
(p s = 0, p s B = 0) which requires simultaneous death of all predators and preys, is not 
dynamically accessible. Thus the 4SPP model can undergo an absorbing state phase 
transition by tuning the different rates with p s B as the order parameter. Our aim here 
is to study this critical behavior in details. 

We have used the standard Monte Carlo methods to study the critical behavior of 
this model. From an initial arbitrary configuration, where each branch A and B are filled 
by arbitrary number of respective particles, the system is allowed to evolve according to 
the random sequential update following the Monte Carlo dynamics of 4SPP model. We 
have studied the system with different values of the rates p, q and r with system size 
L = 10 3 . For illustration, we have fixed two of the reaction rates, say q = 0.02, r = 0.9, 
and vary p as the the control parameter in the following simulations. 

For p < p c the average B density p s B decreases continuously until the system reaches 
a state with no predators (p s = 0) and ultimately the whole prey lattice branch is filled 
with preys. Once reached in this absorbing state (p s = l,p s B = 0), the system remains 
there forever. While, for p > p c the average density of B particles p s B saturates to a 
nonzero value along with p s A which also takes a value smaller than unity. In Fig. [2(a) 
we have plotted p A and p s B as a function of p for a system size L = 10 3 . As expected, 
the density of the preys p s A never vanish and they proliferates in the whole lattice (with 
p s = 1) in the absorbing state where predators are absent. Thus the system undergoes 
an absorbing state phase transition as the birth rate of preys p crosses a critical threshold 
p c = 0.148(4). 

For p > p c , the order parameter p s B shows power law behavior with the distance 
from criticality, 

P S B ~(P-Pcf, (1) 

when p approaches p c . This Eq. ([T|) can be used to estimate p c and (3. As shown in Fig. 
[3(a), p s B versus (p — p c ) is linear in log scale, for the correct choice of p c = 0.148(4); 

| It has been predicted earlier [22] that, in absence of site restriction, the predator and the prey system 
in ID always remain in the co-existing phase. 
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Figure 2. (a) Average steady state density p s A ,p s B is plotted against the control 
parameter p for fixed value of q = 0.02, r = 0.9 and L — 10 3 , critical point is indicated 
at p c = 0.148(4) showing absorbing to active phase transition. Both densities arc again 
plotted against q keepings = 0.55, r — 0.9 fixed; critical point is q c = 0.095(6) showing 
active to absorbing phase transition (see inset), (b) Phase diagram in the (p,q) plane 
for r = 0.9, showing active and absorbing phases. 



the corresponding slope (3 = 0.367(7) gives the estimated value of the order parameter 
exponent. 

One can obtain few other critical exponents from the decay of the order parameter 
p B (t) from an initial configuration with large number of predators. Clearly, in the 
vicinity of critical point p B is a function of time t and the temporal correlation length 
£ M which vary as ~ \p — p c \~ u "- Again, after an initial decay p B (t) ~ t~ a the predator 
density in the active phase approaches the steady state value p s B in the t — > oo limit. So 
p B must scale as, 

p B (t,p) = t-«F(t\p-p c \ v "). (2) 

Thus, one expects that p B (t) for different values of p (shown in the Fig. E^b)) collapsed 
into a single scaling function J 7 , when p B t a is plotted against t\p—p c \ u ". The main figure 




Figure 3. (a) Average predator density (p s B ) is plotted with fixed value of q = 
0.02,r = 0.9; with different choices of p c = 0.153,0.151,0.148,0.146,0.144 shown 
from top to bottom. The correct choice of p c = 0.148(4) and corresponding slope 
in logarithmic plot gives critical exponent /3 = 0.367(7). (b) Time evolution of p B 
with q = 0.02,r = 0.9. Below the critical point for p = 0.142,0.144,0.146 density 
p B eventually extincts and above the critical point p = 0.150,0.154,0.160 density p B 
saturates. At the critical point p = 0.148(4), p B (t) ~ t~ a gives the critical exponent 
a = 0.194(4) and the data collapse gives v,, = 1.8(1). 
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here shows the data collapse when we choose a = 0.194(4), and u tl = 1.8(1). Since at 
the critical point p B {t,p c ) = t _a J r (0), one can obtain both p c and a directly from the 
log scale plot of p B versus t which is linear (as shown in the inset of Fig. EJ^b)). The 
resulting p c and a are consistent with those obtained from the data collapse. Again, in 
the upper critical regime, p s B vanishes as \p — p c \, in the t — > oo limit. This can happen 
only when the off-critical scaling function jF{x) ~ x^l Vn \ thus 

a = p/u„ 

Since all three exponents 0, a and z/„ are calculated independently, one can check if the 
above scaling relation holds. In this case it holds to a great accuracy for the values of 
(3, a and u n calculated here. 

Now we turn our attention to the finite size scaling of p B at the critical point. 
Again, the system of length L with a high density of predators p B (t,L) decays as t~ a , 
indicating a scaling form 

p B (t,L)=t- a g(t/L z ), (3) 

where z is the dynamic critical exponent. Thus, one expects p B for different values of 
L to be collapsed to a single function when plotted against t/L z . This is described 
in Fig. H^a). The inset there shows variation of p B (t) for different system size 
L = 50,100,200,300,400,500, which were made to collapse to a single function using 
a = 0.194(4) and z = 1.52(0). From the scaling relation 

z = f„/fj_, 

one expects that u± = 1.19. This can be verified from the modified scaling relation Eq. 
(J3J). Since za = f3/v±, we have 

p B (t,L)=L-^g(t/L z ), (4) 

where G(x) = xQ{x). In Fig. H](b) we have plotted p B L^I Vl - as a function of tjU 
and found that the data for system size L = 50,100,200,300,400, and 500 could be 




Figure 4. (a) Finite size scaling function p B t a plotted against the scaled variable 
tL- z for L = 50, 100,200,300,400,500 (bottom to top in figures) with p = 0.148(4), 
estimates z = 1.52(0). (b) The same data could be collapsed according to Eq. (J4j) by 
choosing /3/Vj_ = 0.30, which gives an estimate v± = 1-2(2). 
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collapsed into a single curve by choosing (3/v± = 0.30. This gives us u± = 1.2(2), which 
is consistent with the earlier estimation. 

The critical exponents of the model are summarized in Table-H] along with the 
the critical exponents of directed percolation universality class. Clearly the exponents, 
particularly /3 and a, are very different from those of DP-class, indicating that the 4SPP 
model belongs to a different universality class. We have studied the model by varying 
other rates also. 








a 




z 




DP 


0.276 


0.159 


1.733 


1.580 


1.096 


4SPP 


0.367(7) 


0.194(4) 


1.8(1) 


1.52(0) 


1.2(2) 



Table 1. Comparison of the critical exponents between the 4SPP model with the DP 
universality class. 

For example, one can take q as the control parameter, keeping p and r fixed. The 
inset of Fig. E(a) shows the variation of p s A and p s B as a function of q for fixed p = 0.55 
and r = 0.9. Evidently, the order parameter p s B vanishes continuously as q crosses the 
threshold value q c = 0.095(6). The critical exponents obtained in this case was found 
to be consistent with Table-HJ Thus, for any fixed values of r, one expects a line of 
criticality in the (p, g)-plane, which is shown in Fig. Mjo). 

That, the critical behavior of the absorbing phase transition observed in the 4SPP 
model is different from DP, can be visualized from the growth of clusters. The space- 
time diagram, starting from an arbitrary initial configuration is shown in Fig. E£a), 
where the occupancy of species A and B are represented separately in the upper and 
lower part respectively. The evolution of clusters are visibly different from that of 1 + 1 
dimensional directed percolation model. Since the prey species can grow independently, 
the space is always filled locally by preys where predators are absent. Again, regular 
striped structures appear in these figures as the species do not diffuse. A natural question 
would be whether diffusion can drive the system to have an absorbing state transition 
belongs to DP universality class or not. 

In the following we introduce diffusion of both the species explicitly in this model. 
Along with the usual dynamics of the 4SPP model described earlier, both A and B 
particles are now allowed to move to the neighboring available vacant space in their 
respective branches. In this case, it turns out that the clusters evolve more like the 
DP model (see the space-time diagram in Fig. EJ^b)). However, the detailed study of 
the absorbing phase transition reveals that the critical exponents are same as given in 
Table-E 

In presence of explicit diffusion, we choose to study the system with fixed rates 
q = 0.2,r = 0.9. Monte-carlo simulations show that the predator density p s B vanishes 
continuously as p is decreases below a critical threshold p c = 0.323(5). As described in 
Fig. |6]^a), near the critical point, p s B ~ {p — Pc) 13 with /3 = 0.370(9). Again, starting from 
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Figure 5. Time evolution of the 4SPP model for 10 3 time steps with prey A (green) 
and predator B (blue) with system size L = 10 3 for the regime p < p c and at p = p c : (a) 
with asymmetric rules (no diffusion), (b) with asymmetric rules and explicit diffusion. 



a large number of predators, the density p B (t) decays to its stationary value p s B which is 
nonzero only in the upper-critical region p > p c . We find that p B {t) for different values 
of p could be merged to an unique scaling function which satisfy Eq. (T5]) by choosing 
a = 0.190(5) and z/„ = 1.7(5). This data collapse is described in Fig. Mjo)- The critical 
exponents /3, a and v n are more or less consistent with those listed in Table-[U It is not 
surprising that addition of explicit diffusion did not alter the universal behavior. In fact, 
though slow, effective diffusion of predators was already occurring in the 4SPP model 
through the rates q and r. 




Figure 6. (a) Average density of predators (p s ) with q = 0.2, r = 0.9 with 
explicit diffusion: absorbing to active state phase transition at p = 0.323(5) and 
the critical exponent fi = 0.370(9) (see inset), (b) Above the critical point for 
p = 0.330,0.335,0.345 the density p s B saturates and below the critical point for p = 
0.305,0.310,0.315 density p s B eventually extincts. At the critical point p = 0.323(5), 
p B {t) ~ t~ a gives the critical exponent a = 0.190(5) and v n = 1.7(5). 



Some comments are in order here. The fact that the absorbing state phase transition 
in 4SPP model is different from that of DP can be seen from the dynamical mean-field 
analysis (ignoring spatial density correlations). The mean-field densities p % A = (s l A ) and 
P % B = ^ n continuum limit, evolve as, 

8j t = PP^ 1 - Pa) - rp A P B ~ v A 9 -t + B A % (5) 
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a J t= -* lA+TI , lP ,- v °g + D*g, (6) 

where va,b denote the velocities of respective species which appear due to asymmetric 
dynamics, and Da,b are the coefficients of diffusion. The first term in Eq. (jSJ) captures 
the growth of a prey restricted by hardcore interactions and the second term there 
corresponds to the interaction between two species. Again, the first term in Eq. 
(J6j) represents the simultaneous death of two predators. Evidently, these mean-field 
equations have two fixed points: the unstable one (p* = 0, p* g = 0) and the stable 
one (p* = 2^+^' p *b = )■ However, when p B -> faster than p A , one can get 

another fixed point (p* = l,p* = 0) from Eq. (0). This unusual absorbing state 
(jf — l,p s B — 0), as discussed earlier, raises a possibility that corresponding absorbing 
state phase transition can be different from the usual APT to (p* = 0,p* = 0). 
Clearly, for p A = 1, the mean- field equation for the predator density ps (from Eq. 
(E])) is identical to that of the DP. Thus, in higher dimension (larger than the critical 
dimension) one expects that the 4SPP model results in the same mean-field critical 
exponents {fi MF = 1 = p Mi7 ', z MF = 2) as that of DP. Another possible reason for the 
new universality class is the asymmetric dynamics, which generates density dependent 
velocity terms va,b- Note that the absorbing state phase transition in the asymmetric 
contact process [23] belongs to the DP-class, whereas asymmetric updating is a relevant 
perturbation to models with extremal dynamics [23] . Detailed study of the 4SPP model 
with symmetric dynamics will be reported elsewhere. 

In summary, we introduce a two species (predator B and prey A) model in one 
dimension where each lattice site is either vacant or occupied by, a single predator, 
a single prey or both. More than one predator or prey are not allowed at any site. 
The preys are allowed to grow independent of the predators, whereas two predators, 
if present at the neighboring sites, die simultaneously. The species interact through a 
dynamics where the predator produces an offspring by eating a prey from its neighbor. 
When the predator density p s B = 0, even a single prey can invade the whole lattice 
by its independent birth process. Thus, the system has an unusual absorbing state 
(p s A = l,p s B = 0); the other absorbing state (p s = 0,p s B = 0) is not dynamically 
accessible. Using dynamical Monte-Carlo simulation we show that the system shows 
an absorbing state phase transition, as the birth rate of prey p is increased beyond 
a critical value p c keeping death rate of predator q fixed. For a fixed r, the line of 
criticality (p c as a function of q) is governed by a universality class different from the 
most generic one, namely directed percolation. This critical behavior is found to be 
robust against addition of explicit diffusion. Note that the dynamical rules of the model 
is different from other four species predator prey models studied earlier [22] in a way 
that the predators in 4SPP model can not eat the prey at the same site. We believe 
that, non-equilibrium phase transition to a unusual and unique absorbing state, may 
result in different universality class. Further study in this direction could clarify this 
issue. 
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